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Abstract 

The arguments leading to the introduction of the massive Nonsymmetric 
Gravitational action are reviewed [Q,^] , leading to an action that gives asymp- 
totically well-behaved perturbations on GR backgrounds. Through the anal- 
ysis of spherically symmetric perturbations about GR (Schwarzschild) and 
NGT (Wyman-type) static backgrounds, it is shown that spherically sym- 
metric systems are not guaranteed to be static, and hence Birkhoff 's theorem 
is not valid in NGT. This implies that in general one must consider time 
dependent exteriors when looking at spherically symmetric systems in NGT. 
For the surviving monopole mode considered here there is no energy flux as 
it is short ranged by construction. Further work on the spherically symmet- 
ric case will be motivated through a discussion of the possibility that there 
remain additional modes that do not show up in weak field situations, but 
nonetheless exist in the full theory and may again result in bad global asymp- 
totics. A presentation of the action and field equations in a general frame is 
given in the course of the paper, providing an alternative approach to dealing 
with the algebraic complications inherent in NGT, as well as offering a more 
general framework for discussing the physics of the antisymmetric sector. 
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I. INTRODUCTION 



In General Relativity (GR), just as in Maxwell's electrodynamics, one finds that given a 
spherically symmetric system, there are no dynamical degrees of freedom in the theory This 
is Birkhoff 's theorem, and implies that a time dependent source will not excite modes in the 
gravitational system, so that outside this source the system must be physically equivalent 
to the Schwarzschild solution (the non-trivial static spherically symmetric spacetime). It 
has been further established |3],|J that the solution is stable when perturbed, so that small 
deviations from spherical symmetry do not alter the large scale features of the spacetime, 
and systems that are only approximately spherically symmetric are therefore still very well 
modeled by the Schwarzschild solution. This establishes that the phenomenology of the 
Schwarzschild solution is physically relevant. 

To see that Birkhoff's theorem is not a generic feature of physical theories one need 
look no further than a scalar field. However the example that is important for this work is 
the (massive) Kalb-Ramond action || which (as will be demonstrated in Section [TTT| ) has a 
single mode that is time dependent in general. As will be shown in Section ||, the massive 
Nonsymmetric Gravitational Theory (mNGT as opposed to the older versions of the 
theory, referred to as NGT, or massless NGT) becomes identically a massive Kalb-Ramond 
field with an additional curvature coupling term when considered as a perturbation about a 
Ricci-flat GR background, so the result that NGT has a monopole mode is not surprising. 
The mode considered here is short ranged, so that far enough away from the source one 
finds that the solution will be dominated by Schwarzschild behavior, and there is no energy 
flux. However, after demonstrating that spherically symmetric fields in the skew sector are 
not static in general, it will be shown in Section [VI| that the symmetric sector will also no 
longer be static, through examination of a similar perturbation about the mNGT background 
discussed in Section [V]. This means that no static solutions can be considered rigorously as 
an exterior solution unless the solution is globally static (i.e. the interior is static as well). 

The results in this paper are obtained through an examination of linearized perturbations, 
although the conclusions must hold in general as Birkhoff's theorem would imply that these 
fields must be static as well. What cannot be examined in this fashion is whether in the 
full nonlinear theory more modes become excited. In particular, one will see in the case of a 
perturbation about a GR background (in Section |IJ) that there are three propagating modes, 
even though the absence of gauge invariant kinetic terms in the full theory would suggest 
that all six modes in the antisymmetric metric could be independent degrees of freedom 
(as yet, the number of degrees of freedom in mNGT has not been rigorously established). 
Although this issue is not addressed directly in this work, the ability to recast the theory in 
a general basis given in Section [TV] sets the stage for a complete analysis of the spherically 
symmetric system in NGT. Given that there is an additional mode in the general spherically 
symmetric system, how the fields may or may not approach an asymptotically flat spacetime 
should then be addressed, and also whether evolution singularities of the type discussed in 
M are encountered. 
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II. MASSIVE NGT 



The original version of NGT [|7| grew out of a re-interpretation of the Einstein-Straus 
HU unified field theory as a purely gravitational system. The antisymmetric part of the 
metric and connection operationally produce different modes of parallel transport and index 
contraction [^|.|10|.|TTJ , where the algebra is consistent with an enlargement of the tangent 



vector space to its hyperbolic complex extension [|i2|4T4]|. It is important to note that the 
action cannot support the additional Bianchi identities and gauge invariance related to the 
extension of the tangent bundle, simply because the base manifold is locally diffeomorphic 
to TZ 4 , and the variational principle is based on an integration over this real manifold. Any 
change of gauge that mixes real and hyperbolic complex covectors will cause the volume 
element to pick up a hyperbolic complex piece, and the action will no longer be real. 

The hyperbolic complex structure is unnecessary for an operational discussion of the 
theory (although it may be relevant for a more fundamental discussion of its physical inter- 
pretation), and in this paper all quantities will be considered real, allowing antisymmetric 
contributions to the metric and connection coefficients. The dynamics of the theory will be 
determined from the first order action (G = c = 1): 

S = J ^{-g^it^T] - grdfrW v] + 1T„ + \a%^W,W v + im 2 g^ [H } + S M , (1) 
where: 

8S M 



TV, (2) 



is the matter stress energy tensor that acts as a source in the gravitational field equations. 
The Ricci-like tensor for NGT appearing in (|I|) is written as: 

and a mass term for g< i has been included along with a term quadratic in W as new features 
of the action. As will become clear shortly, a may be fixed uniquely by requiring good 
asymptotic behavior of perturbations about GR backgrounds. The new parameter in the 
massive action (m) is an inverse length scale that must be constrained by experiment, and / 
is a Lagrange multiplier employed to enforce the vanishing of the trace of the antisymmetric 
part of the connection coefficients: = TfL^. 

The field equations related to the metric compatibility conditions are derived through 
the variations: 



6S_ 

SS 

5S 



r„ = 0, (4a) 
d v gl v>l] + ag^W u = 0, (4b) 
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Contracting on either index of ( f4c|) and solving for the (anti-) symmetric parts of the diver- 
gence of the densitized inverse metric results in the determination of the Lagrange multiplier 
using (|4"T5|): 1°" = ^g^^W^. This also allows one to simplify the Kronecker-5 terms, and 
determine the compatibility conditions in undensitized form as: 

9 7 g^ - gfmK-y - 9auT^ = |a(^[ 7 p a ]i/ + l^/w^a])^ ^^/?, (5) 

where the inverse of the metric has been defined by g lia g ai ' = 5 V , which has been used in 
order to rewrite the compatibility conditions in terms of the components of the metric g^ v . 

The remaining field equations derived from the variation of the action with respect to 
g^ u may be written as: 

K, u := R™ + d^W v] - \aWpW v - \m 2 (g [H - g a tf v0 gW + \g„»g [aP] g [aP] ) = T„ v - \g vtl T, 

(6) 

where T = g^T^, and the tensor 1Z has been introduced in order to simplify the discussion 
of the field equations. One may translate the conventions used here to those in [|TJ by taking 
W — > —^W,T^ U — > SnT^, a — > — |cr, and adjusting the definitions of the inverse metric: 

— ¥ g u ^. To see the equivalence of the action, one further needs to rewrite T in terms of 
the unconstrained W connection, and drop the contribution from the Lagrange multiplier I. 

The action for massless NGT is given by ([D with m = a = 0. As will be demonstrated, 
the new terms have been introduced in order to make all skew modes short ranged when 
considering perturbations about GR backgrounds. One performs this expansion about a 
symmetric, Ricci-flat background, where one assumes that all background curvatures fall 
off at worst as 1/r as r — * oo. This allows one to talk sensibly of energy-momentum 
and decompose fields via a spin projection, so that higher order poles and negative energy 
(ghost) modes may be identified, as well as avoiding the full nonsymmetric structure of a 
more general background that would make the analysis far more complicated. One considers 
a perturbation of all quantities about a symmetric GR background as in ||15|| : 

g^v * 9(fj-u) hfj,v> 

1%,- {£} + 1% (7) 

where W, I and T are are considered to be first order in the perturbation (as the background 
is assumed to be Ricci-flat, there one has T^ u = 0). As usual, indices will be 'raised' and 
'lowered' by the symmetric background metric, and the covariant derivative V is associated 
with the background Christoffel symbols { } determined from the background metric in the 
usual manner. Corrections to the background curvatures, field equations and Lagrangian at 
each order in the perturbation will be indicated by a superscript as: °R, l R ■ ■ -. 

The first order correction to the compatibility equation (||) can be solved explicitly for 
7 to yield: 

1% = lg aP (Vvhp„ + VA/J - V^V) + 1<*6%W V] , (8) 
and 7 M = 7^ is seen to vanish by the linearization of the skew divergence equation: 

_ V^ [H = aF. (9) 
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In massless NGT ||10|| , one had (§,^D with a = 0, and hence there was no relation between 
the metric degrees of freedom and those of W. The skew part of the linearization of equation 
(^) with m = as well as a = became: 

%h = X R[M + d [fl W u] = V a7 f H + d^W v] = T [H . (10) 

The symmetric contribution, are the equations for a metric perturbation in GR and will 
be ignored in the remainder of this section. Using (|9]) with a = 0, the assumption that the 
background is Ricci-flat ( i? M ^ = 0), and the commutation relation (for an arbitrary tensor 
B): 

V^V/jj-B^ = -^{B^R ^ + B^ R"^), (11) 

(|TCi|) simplifies to: 

V Q V a /i [H - 2d [lx W u] - A R a /uh[ a( s] = V a F^ a - 2d^W v] - 8°iT/> M = -2T [H . (12) 

The second form is given in terms of the curl of the skew metric {F ltiv = d^h^ u + <9„/i 7M + 
d^h^), in order to more easily demonstrate the result found previously by Damour, Deser, 
and McCarthy [|T^,|T^] (using the fact that V 7 V^F Atiy7 = about a Ricci-flat background), 
that these antisymmetric perturbations will in general have bad asymptotic behavior. Al- 
though written as if the skew metric were a gauge field, the presence of the curvature coupling 
term implies that the associated gauge invariance is not present [ I9| , and one is not allowed 



to make any choice of gauge in order to simplify this sector. One proceeds by taking the 
divergence of (|12D and choosing the gauge V a W a = (the theory had a U(l) invariance as 
W only appeared in the action in a curl) to find: 

V Q V Q ^ - SVT^VlV/d = -2V^T [H . (13) 

Notice that the background curvature acts as a source here, so that even if one postulates 
that the matter source is conserved, this curvature coupling (and in general other nonlinear 
terms) will still exist as a source, causing W to propagate and have asymptotic behavior 
consistent with a massless field (~ 1/r along the forward light cone). Using this asymptotic 
behavior to determine h from ( |I2"D results in a source with ~ 1/r behavior, causing the field 
h not to fall off as r — >• oo along the forward light cone. This analysis is correct since one 
has assumed that the background curvature falls off fast enough, and hence the potential 
term in ([13|) can be treated as a source, without changing the asymptotic behavior of the 
fields. It must be stressed that one is assuming that the background and the radiative fields 
fall off at least as fast as ~ 1/r, and what has actually been derived here is a contradiction 
of this, since h is driven to a constant and can no longer be considered as a perturbative 
mode. 

Any analysis of this sort also supposes that a solution of the linearized field equations 
does in fact correspond to an exact solution of the full nonlinear field equations. This 
is the case in GR [|T6i^0| at least for source free equations, but no such result exists yet 
for any of the models considered here. It is possible to take the stance that NGT is not 
linearization-stable, so that this sort of analysis necessarily produces spurious results that 
do not correspond to global solutions, but then one is denying the ability to do any sort of 



5 



perturbative analysis without the existence of an exact solution to back it up. Due to the 
scarcity of solutions, and the apparent existence of weak-field perturbative situations, this 
would seem an unreasonable position to adopt. 

This result is not confined to curved backgrounds, and in fact the analysis about 
Minkowski space will serve to explicitly demonstrate the higher order pole leading to bad 
fall-off. Since this curvature coupling, and any nonlinear effects in general, will act as a 
nonconserved source term in the skew sector, the linearization that correctly represents the 
full nonlinear field equations in the asymptotic region will have a nonconserved source term. 
This is no more than the observation that once again, the full NGT action does not possess 
any form of additional gauge invariance in the skew sector. A gauge field coupled to a source 
(or matter) in a non-gauge invariant manner may have drastically different behavior than the 
empty space and apparently gauge invariant field equations, if indeed the action is consistent 
at all. Any analysis that attempts to determine the propagator or asymptotic behavior of 
the field must take the form of the source (or coupling to other fields) into account. 

A trivial example of this is given by considering the Maxwell action. Coupling the 
usual gauge invariant kinetic terms to a nonconserved source gives an inconsistent set of 
field equations, and adding some sort of gauge fixing term will give a consistent set of 
equations, but the scalar ghost mode will be excited and depending on the gauge there may 
be higher order poles in the solution. The linearized field equations considered outside the 
source resemble those of the gauge invariant theory in a particular gauge, but treating them 
as such will not give asymptotic behavior that follows from coupling to the nonconserved 
source. The situation in NGT is more akin to enforcing the gauge condition in the action 
through the use of an auxiliary field as: bd^A 11 [^1],^]. Source conservation and absence of 
ghosts relies on whether or not the scalar Lagrange multiplier field b has a source or not in 
the wave equation that determines it, and is thus a global question. Given that the source 
for A is not conserved, then b propagates and there are higher order poles in the solution 
for A, leading to fields that do not fall off as r — ► oo along the forward light cone. 



Considering the field equations for massless NGT linearized from (12) about Minkowski 
space: 

Dh [H - 2d^W u] = -2T [H , (14a) 
d u h [ ^ ] = 0, (14b) 

(d u T[ UfJ i ^ 0), one may take a divergence (to find a wave equation for W) or a curl (to remove 
W) of the first of these, resulting in the unique consistent solution: 



h^] = —2D 1 



T\jjj,] + 2D 1 d a d^T{ au }} 
-lav, 



w M = -2n- 1 a"T M . (is) 

The presence of the higher order pole (and consequent bad fall-off) is now obvious from 
the presence of the D~ 2 term in the Green function solution, and is no more than the 



result of vanNieuwenhuizen |23[] who showed that the only healthy quadratic actions built of 
antisymmetric tensor fields are the so-called Kalb-Ramond (massless or massive) actions. 
One also notes that there are 5 modes here: 3 in h, since 3 are determined algebraically 



by the second equation in (14]), and 2 in W due to the previously mentioned U(l) gauge 



invariance p4J . If it is assumed that the (matter) source is conserved, the higher order poles 



6 



are removed at linear order, but show up in the second order correction to the fields, again 
causing a breakdown of the perturbative analysis. 

This analysis correctly represents the asymptotic behavior of the fields (W, h), and is 
equivalent to equation (18) of p5| , where the higher order pole resides in the projection 
operator: -P(l + ). One also sees the true propagating nature of W, and this is borne out by 



the analysis in pq , |27 l where there are five degrees of freedom evolving from each Cauchy 
surface, the extra two of which are associated with the field W. That a Lagrange multiplier 
is propagating merely signifies that it is a determined multiplier, with its evolution derived 
from the field equations jr5| and not freely fixable as was done in PB"| , |25| and in the next to 
last section of [^J where ad hoc constraints were imposed on the linearized theory in order 



to obtain the dynamics of a Kalb-Ramond theory. That these constraints cannot exist is 
clear from the lack of gauge invariance in the full NGT action. 

The result of vanNieuwenhuizen does however motivate a potential solution to this prob- 
lem, since the massive Kalb-Ramond theory does not require a conserved current and yet 
has no ghost modes, higher order poles or tachyons. The additional terms in the action for 
mNGT (|1|) are introduced in order to allow the linearized field equations of NGT to take on 
this form in the antisymmetric sector. These two terms play slightly different roles: the W 2 
term causes W to be determined in terms of metric functions directly (a is fixed in order 
to find the correct form of the kinetic energy terms), and the mass term for gn makes the 
skew sector short-ranged, and ensures that the linearized field equations remain consistent 
when expanding about a flat background. 

Thus mNGT should have a linearization about Minkowski space of the form: 

d a F ilva + m 2 h [liu] = J [H . (16) 
The solution to fllTf) can be found by taking a divergence and substituting back in to find: 



h[pa/] = 



Phw = ^J M . (17) 

The higher order poles have disappeared, and it can be shown that the linearized Hamil- 
tonian is weakly positive definite and that ghost modes are removed through the algebraic 
conditions that couple them locally to the source in fllTl) . About a more general background 
one can allow a curvature coupling term, since it will not affect the behavior of the fields 
asymptotically once the background is assumed to fall off appropriately. Choosing the theory 
that results in this behavior in the linearized theory will fix a uniquely. 

Returning now to the field equations of mNGT expanded about a GR background fol- 
lowing from one finds: 

i^iiv = R^v + ^V[ M V a /i[ ai ,]] — \m 2 h[^ = T^ u — lg Ufl T, (18) 
where the first order correction to the 'Ricci' tensor is given by: 

%„ = V a7 * -V^)- (19) 



Again ignoring the symmetric GR perturbations, the antisymmetric part of (18) is: 
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= -\{V a F^ va + m 2 h [H ) - 2V a V [fl h [au]] + I[l + |a]V [M V"/i M] = T [H . (20) 

Requiring that this reduce to the massive Kalb-Ramond field equations (|I6|) determines the 
(previously arbitrary) coupling: a = 3/4. The last two terms can be reduced to a curvature 
term to give: 

V a F^ a + m 2 /i [H - 4°ir t/ f v] h [afii = -2T [H , (21) 

so the skew sector perturbations are well-behaved when perturbing about any asymptotically 
flat GR background. 

Expanding the action (|1|) to second order (ignoring surface terms) gives: 

2 C = - 2 R - \h}R + h^% u + h^d [fl W u] + laW^W^ + /T M - \m 2 h^h [tiv] (22) 

and once compatibility is imposed, followed by the removal of W, this becomes: 

2 C = ±F^F^ - lm 2 h^h [H - VhMV„h bll] - + ±)V„/^W% 7] . (23) 

Choosing a = 3/4 results in kinetic terms identical to those of Kalb-Ramond theory on a 
GR background, giving the skew sector action: 

2 C = ±F^F^ - \m 2 h^h [H - h^h^°R a , Pu , (24) 

which reproduces the linearized field equations (|2T|). Thus the massive NGT action will be 
(H) with a = 3/4 [|lj , giving the action (|4]) for perturbations about a GR background and 
guaranteeing good asymptotic behavior for these fields. 

Although it has been established that the perturbation equations about a GR background 
are a consistent system resulting in good fall-off for the skew sector, it is not clear whether 
an asymptotic perturbation actually corresponds to a global solution (linearization stabil- 
ity). The (seemingly contrived) asymptotic limit of gauge invariant kinetic terms cannot be 
reflected in the full action, since there is no room for the additional gauge invariance in the- 
ories constructed from antisymmetric fields in this manner. This means that in general that 
one expects more (perhaps all 6) degrees of freedom in the skew sector evolving as degrees 
of freedom in a Cauchy analysis, whereas in any spacetime that has an asymptotically flat 
region only three will survive. This situation could be similar to that found in 0, where 
vector fields were seen to increase their degrees of freedom when gravitational effects are 
taken into account. In order to obtain an asymptotically flat spacetime (with the reduced 
degrees of freedom of the vector fields) from physically reasonable initial data, the evolution 
equations were seen to have to encounter singularities. This is generally considered to be a 
sign of instability, and certainly not a desirable feature in any theory. Perturbations about 
NGT backgrounds should also be considered, since the physically interesting NGT solutions 
are most likely those that are not 'close' to a GR solution |3T|| . It is hoped that a more com- 
plete analysis of the general spherically symmetric system should be able to say something 
about this issue, since it will certainly tell one how many degrees of freedom survive and 
how they couple to external fields, and hopefully something about how the system may or 
may not approach an asymptotically flat spacetime. 
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It is also true that the form of the action (||) is far from unique. In particular, one could 
replace the W 2 term with some combination of W 2 and :y=g( M z,)d Q [g^ a ']d 7 [g^], giving 
the same perturbation equations (|2~lf) , and resulting in an arbitrary coupling constant in 
the action. Further, since there is nothing preventing one from adding T[] terms (they are 
tensors) or even infinite strings of terms of the form: 5 , ( Q/ g)5'^ 7 ' ) ■ • ■ or gi a p]g^ • • • (which 
conveniently disappear in the asymptotic expansion), there is clearly an infinite number of 
actions that do this. These examples seem extremely unnatural and will not be considered 
further here, although the results of this paper would not change significantly for these more 
general actions. 



III. SPHERICALLY SYMMETRIC PERTURBATION OF THE 
SCHWARZSCHILD SOLUTION IN A COORDINATE BASIS 

The absence of a Birkhoff theorem may be derived from the perturbation equations fl2~T|) 
developed in Section |I[ In general the spherically symmetric fields in the skew sector will not 
be static, although the symmetric sector will remain static in the perturbation about the GR 
solution considered here. The background metric is Schwarzschild with (coordinate basis) 
metric written as: g = diag(A(r), —1/A(r), —r 2 , — r 2 sin 2 (#)), where A(r) = 1 — 2M s /r and 
M s is the Schwarzschild mass parameter. The perturbation considered will be one that is 
spherically symmetric but not necessarily static. A killing vector analysis yields the general 
form of the spherically symmetric perturbation: 



h 00 (t,r) h ( p 1 )(t,r) + h[ i](t,r) 

/i(oi)(t,r) -/i[ i](t,r) h n (t,r) 

h 22 (t,r) h m (t,r)em(6) 

-h [23] (t,r)sm(9) h 22 (t, r) sin 2 (#) 



(25) 

Making a change of coordinates of the background geometry is equivalent to making a change 
of gauge on the perturbation: 5h = £ e [g], where e is the spherically symmetric vector gauge 
parameter generating diffeomorphisms between spherically symmetric spacetimes. This al- 
lows one to simplify the form of the perturbation by a suitable choice of gauge. Choosing 
the gauge parameter as: 




d t [h 



22 



2rA{rf 



dr, (26a) 

(26b) 



removes the 9 — 9, <fi — <fi, and symmetric t — r perturbations altogether, and a remaining 
gauge transformation e° = e(t) allows one to remove an arbitrary function from the t — t 
component of the form: 5h 00 = 2A(r)e(t). 

The field equations will be written without the source terms for simplicity although it 
is straightforward to include them and relate the constants of integration to properties of 
the source. First reviewing how the symmetric (in this case identically GR) perturbations 
become static, it is simplest to begin with the field equation: ^(oi) = 0, which implies: 
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d t [h u (t,r)] = 0, 



(27) 



immediately showing that hn must be static. By considering ^22 = 0, it is determined to 
be: 

M*,0 = -^, (28) 

where the integration constant has been combined with M s and interpreted as a perturbation 
of the Schwarzschild mass parameter: 8M S . Then one considers: g 00 ^oo — 9 11 ^11 = 0, 
leading to: 

h 00 (t,r) = B(t)A-^, (29) 

r 

also giving a contribution arising from the perturbed mass parameter, as well as an arbitrary 
function of time as an integration constant, removable by the remaining choice of gauge noted 
above with: e{t) = —B(t)/2. Thus one has that the symmetric perturbations are static and 
interpretable as being due to a small change in the total energy of the system: SM S . 
In the skew sector, the t — r field equation gives: 

1 ^ [01] = |(2A // -m 2 )V] = 0, (30) 

from which one must conclude that h[oi\ vanishes outside the source. (Primes will denote 
the derivative of a function of one variable where convenient.) In massless Kalb-Ramond 
theory, this is the surviving spherically symmetric ghost mode which in that case is pure 
gauge. When a mass term is added, although these modes are now no longer pure gauge, 
they do not propagate since they are locally coupled to the source. It is these modes that one 
eventually must worry about, since in the full theory they may play a nontrivial dynamical 
role. In the 9 — sector: 

1 ( 1 24 4 1 

^[23] = - J [^[h m \ - Ad 2 r [h m ] - (A' - —)d r [h m ] + (-(1 - A) + m 2 )h [2 ^sm(9) = 0. 

(31) 

In order to derive the asymptotic form of this perturbation, it is convenient to define h\23\ = 
rf(t, r), leading to: 

-d 2 t [f] - Adl\f] - A'd r [f] + (— + ^2- + m 2 )f = 0. (32) 



Introducing the coordinate: 



(h r + 2M,ln(-£--l), (33) 



Mr) ' 2.1/.. 

one obtains (after multiplying by A) the partial differential equation for / in normal form: 



9 A 1A' 

dUf] - d 2 Af] + A(m 2 + ™ + 2±)f = (34) 
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where r is considered as a function of r* as are A(r) and d r A(r), and the perturbation 
/ = f(t, r*). In this form it is obvious that (0) is a hyperbolic wave equation, and that the 
field / is therefore nonlocally related to the source. 

Using the fact that 1/r — 1/r* ~ o(l/(r*) 2 ) as r* — > oo, one keeps only the constant 
mass term asymptotically, as all other potential terms will be dominated by it. This leaves 
the massive scalar wave equation to determine the asymptotic form of the perturbation: 

dm-d 2 M+m 2 f~0. (35) 

The static solution of this is easily seen to have the asymptotic form: 

h m (r) ~ F -e- mr * ~ F -e~ mr {^rY 2mM % (36) 
m m 2M S 

where a factor of m has been introduced in order to make the constant Fq dimensionless. 
The general time dependent case may be handled by noting that the retarded and advanced 
Green functions for the massive scalar wave equation plj (x 2 = t 2 — x 2 ): 



jyr<&flAv, x \ = D r e t,advU r \ = — 0(±X°) 



2 m6(x 2 ) i-z. 
5(x 2 ) ^Ji(mVi 2 ) 



2Vx 



(37) 



depends only on (t,r), and r *£) ret ' adv (t ; r *) w ill solve d35|) . The asymptotic behavior of h[2g\ 
is then determined from: 



l[23] 



rr*D ietMv {t,r*). (38) 



Note that the behavior on the light cone is determined from just the massless Green function 
S(x 2 ) and so it would appear that hm] will behave as r as r — > oo along the forward 



light cone. This is misleading, as it can be demonstrated explicitly |33|] that for C°° initial 
data with compact spatial support, a massive Klein-Gordon field is bounded everywhere by: 
< d(l+ | t |) -3 / 2 , for some constant d, and therefore cannot radiate energy. This can 
also be understood by noting that because the field is massive, the effects propagating on 
the light cone must be fields of infinite energy, and given some physically reasonable source 
distribution, these infinite energy modes will not be excited. 

The existence of time dependent solutions thus proves that Birkhoff theorem is not 
valid in mNGT, although the short-ranged nature of the skew sector implies that monopole 
radiation will not exist. The symmetric sector has remained static in this system, but as will 
be shown in Section [VT] , through a perturbation about an approximated mNGT solution, 
this will not be the case in general. The perturbation equations about a mNGT background 
have not been given in covariant form, primarily due to the complication involved (although 
it is possible in principle using a generalization of the inversion of the compatibility equation 
given in p4]| ) . Instead the system may be developed in each case separately, and the analysis 



simplified by considering the field equations in a vierbein frame given in the next Section. 



IV. NONSYMMETRIC THEORIES IN A GENERAL FRAME 

The structure of the compatibility relations and field equations in nonsymmetric theo- 
ries can be formulated in terms of components in a general moving frame (in the sense of 
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global section of the general linear frame bundle GLA4 of all linear frames over M). The 
formalism given here is essentially a more systematic development of the approach in ||35|| . 
and differs slightly from that of Hlavaty |36|] in that the (in general nonsymmetric in a co- 
ordinate basis) connection coefficients have been split up into a connection that is torsion 
free, and another that is purely antisymmetric, instead of defining two types of covariant 
derivative, one associated with the NGT Christoffel symbols, and another that is in general 
non-symmetric and not in general torsion-free. The construction here has the advantage 
of only defining one covariant derivative, and the fact that it is torsion-free implies that 
the antisymmetric components in a general (non-coordinate) frame are related in the stan- 
dard way to the structure constants. In a coordinate basis this is the usual split between 
the symmetric and antisymmetric components, however it is easily generalized to any basis 
by considering the antisymmetric components as a separate antisymmetric tensor, and the 
symmetric components as a torsion-free but generally non-compatible connection. 

This provides a simple way to split the GR and NGT contribution in weak field situ- 
ations, as well as generating computationally simpler systems to solve when inverting the 
compatibility relations. Note that although the formalism is developed for a general basis, 
the specialization to a vierbein basis (the reduction of GLM. to LA4, the Lorentz frame 
bundle consisting of all Lorentz frames above M) which will be utilized in the rest of this 
paper, is accomplished through the choice of the fiber metric as ) — > r] above all points of 
the manifold. This is possible in NGT for the same reason that it is possible in GR: math- 
ematically formulating a physical theory in a diffeomorphism invariant manner will always 
allow the introduction of these general linear frames. The reduction to Lorentz frames is also 
possible as one is assuming that the symmetric part of the metric that one is attempting to 
diagonalize is nondegenerate, allowing the reduction of the frame bundle. This construction 
will be of importance when considering the canonical analysis of NGT, as one would like to 
work in a surface compatible (generally non-coordinate) basis in order to avoid specialization 
to a particular choice of time parameter fixed by the foliation of the manifold, and is easily 
applied to other systems with a nonsymmetric metric and connection plj . 



A. Metric, Compatibility and Curvature 

The compatibility conditions in a coordinate basis (|5p will be written for convenience as: 

where A depends only on the metric or quantities directly derivable from it (and possibly 
other quantities, but for the purposes of this construction it does not depend on the con- 
nection coefficients). Parallel transport (and the related covariant derivative) will then be 
defined using just the symmetric part of the coordinate basis connection, and its action on 
the (coordinate) basis vectors is: 

Ve. [e], = r? Q/3) e 7 , V eQ W = -Vj ap) 9^ (40) 

and the connection is split into a symmetric connection and an antisymmetric tensor: 

r^-r^+A^,. (4i) 
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Thus T will refer from this point onwards to the torsion-free (symmetric in a coordinate 
basis) part of the connection, and A to the remaining tensor contribution. In this way, T is 
a torsion-free (but non-compatible) covariant derivative since: 



T 7 = Q~ 1 



= 0. (42) 



The compatibility equation ([39]) then becomes: 

^e 7 [fi , ]/^i/ = &-y{9ftv] dceu^-yfi = + S«i/A^ ^j/i (43) 

where the basis vectors are just directional derivatives along the coordinates: e 7 [ ] = d 7 []. 

With this definition of the covariant derivative and related connection coefficients, the 
geometric curvature is found as usual from: 

R%» = ^[(Ve M V e „ - V e „V eM - V M )ep) 

= e»[K p ] - e u [T« p ] + TZ p r>„ - T^, (44) 

and defining the two independent contractions: 

K» = R %» = - e ^ + r ^ r «7 - (45a) 

The Ricci tensor will be defined as: 

p _ pi 1 p2 

= e a{^ufi\ ~ 2 Fa/J + e A*P-\, a ]) + ^ZfJ^aj ~ ^a/J^v-y' (46) 

This particular combination is symmetric, and obviously reduces to the GR Ricci tensor 
when the NGT antisymmetric terms vanish. Decomposition of @ into R^ v and another 
that depends on A as: = R^ u + R~, where (A^ = A^ Q ) gives: 

Riu = V eQ [A]£, + V £[m [A]„] + A^Af a . (47) 

A more general basis is introduced at each point on the manifold through e A = Ej^e^, 
where E is locally an element of Gl(A,TZ), and these bases are smoothly joined up to form 
sections of the tangent bundle T(A4) |37| , |38[| . The general basis vectors are then given 
in terms of a coordinate basis through the vierbein-like quantities, which can be used to 
translate tensors from one choice of basis to the other: 

e A = E/e^ E/E/g^ = g AB , etc.. (48) 

(In the usual orthonormal basis, one transforms the symmetric part of the metric to the 
Minkowski space metric, and the £"s provide the isomorphism between coordinate basis 
tensors and locally Lorentzian tensors.) The dual basis of T*(M) is introduced through the 
usual relation: ^[ee] = 5g, and the inverse of the vierbeins is defined through: E^E A V = 5%. 
In this paper capital letters from the beginning of the alphabet: A, B,C, ■ • • will refer to 
components of the object decomposed in the general basis. 

Parallel transport of the basis vectors now defines the generalized connection coefficients: 
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V eA [e] B = r%ec, \/ eA [0f = -T C AB e B . (49) 

The definition of the basis in ( |4*BD implies that it is no longer a coordinate basis in general, and 
hence the directional derivatives no longer necessarily commute, giving rise to the structure 
constants: 

[e A , e B ] = C AB C e c , (50a) 

given by: 

C AB C = E c u (E/d,[E B »] - E/d,[E/]), (50b) 

calculated by noting that e A [ ] — E^d^l ]. This also implies that a torsion- free connection 
will no longer be symmetric, and vanishing torsion now gives: 

Tic = A [V eB e c - V ec e B - [e B , e c ]] = 2Tf BC] - C BC A = 0, (51) 

allowing one to determine the antisymmetric part as usual from the structure constants. 
This is the motivation for splitting up the connection in this way. Given some alternate 
split where V is not torsion free, one would have to distinguish between the effects of the 
general basis on the skew part of the connection coefficients, and that of the NGT effects 
(themselves tensors). 

The compatibility condition (|39"D can now be written as: 

V ec [g) A B = 9adA bc + 9dbAEa ~ &cab, (52) 

where since A and A are tensors, they are just redefined by multiplication by the appropriate 
combination of vierbeins. The symmetric part of this can now be solved for the symmetric 
part of T in terms of the antisymmetric part, the structure constants, and A, to give: 

r<7(AB) = \&-C{AB) — ^A[BC] + ^B[CA] — A D A A DBC + A D B k DC A, (53) 

where the quantities: 

^abc = g(AD)F B c, Aabc = 9(ad)A bc , A A B = S^ AC ^ g[ C B], (54) 

have been defined for convenience, and S is the inverse of the symmetric part of the metric 
defined by: S^ AB ^g^ B c) = S A . Also appearing is the symmetric part of: 

Acab = e B [g CA ] + e A [g BC ] - e c [g AB ] + A° BCA + A° ABC - A° CAB . (55) 

The antisymmetric part of the compatibility conditions can now be recast (using (|53|)) as 
24 algebraic equations for A: 

Acab - A D A A E B A ECD - A D A A E C A EBD + A D B A E A A ECD + A D B A E C A EAD = Q C[AB] , (56) 

where: 

^C[AB] = |(A<7[AB] + A D B A D ( CA) — A D A A D ( BC )) 

+ A D A (T B { CD ] + T C [ BD }) — A D B {T A [ CD ] + T C [ad}) — A D c T D y AB y (57) 
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The method for solving the compatibility conditions is to first determine the auxiliary 
quantities appearing in this relation: (A, Fn,A,Q) in terms of the vierbeins and metric 
quantities, then solve for A through (|56|), determine T( ) from (|53"D, and then use S with 
and rV ) to form T BC and A BC . This may not seem like much of a simplification, but when 
specialized to a Lorentz frame, many of these quantities simplify considerably (as is the case 
in the Wyman sector in Section |V|). 

The curvature tensor ([0]) becomes: 

rA bcd = # A [(V ec V eD - V eD V ec - V [ec>eD ])e jB ] 

= e c[Fr> B ] ~ e £>[rce] + Tdb^ce ~ ^cb^de ~ Ccd E ^ebi (58) 
and the contractions: 

R AB = RC ACB = e c[^BA\ - g b[^Ca] + ^BA^CD ~ ^CA^BD ~ CcB^DA^ ( 59a ) 

r ab = R cab = e A[^Bc] ~ e B[^A.c] - C AB D T DC , (59b) 
combine to give the Ricci tensor: 

p _ pi 1 p2 

J^AB — J^AB ~ 2 AB 

= ec[T C BA ] - e B [Y c CA ] - \e A [Y c BC ] + \e B [T c AC ] 

, pD r C r D r C ri D r C i l/o D r C 
* BA CD ~ 1 CA L BD ~ ^CB J- DA ' 2^AB J- DC 

= e c [T c BA ] - e B [Y c CA ] - \e A [T c BC ] + \e B [Y c AC ] 

+^BA^CD ~ ^CB^DA + ^\AB]^DC- (60) 



In the split Rf B = R AB + R\ B 

V rAl<7__uY7 



R% = V ec [A]2b + V e , A [A]bj + A AD A BC (61) 



as expected. 

Since by construction T is a torsion free connection, and (|58|) is the standard curvature 
tensor constructed from it, one obtains the usual Bianchi identities [^8| on the curvature 



tensor. One should note though that the connection is not compatible, and so the rotation 
coefficients are not antisymmetric. The relevant Ricci tensors are also not constructed in 
the same manner as in GR, so the implications of these identities are somewhat different. 
The first Bianchi identity gives the usual cyclic identity on the last three indices of (p8|). 
and leads to the result: 

R[ab] = \R\b (62) 

when one contracts on any lowered index. (This can also be proven directly using the Jacobi 
identity.) This tells us that the NGT Ricci tensor is symmetric {R\ AB ] = 0) in general, not 
just in a coordinate basis. 

A detailed study of the contractions of the second Bianchi identity (the cyclic covariant 
derivative): 

V ec [R] a BD e + V eD [R] A BEC + V eE [R] A BCD = 0, (63) 



should result in a derivation of the equations of motion for matter fields p9[ , |40"f from the 
field equations. 
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B. The NGT Action and Field Equations in a General Basis 



The translation of the field equations (f|fj,|5L[3) is accomplished through an almost straight- 
forward substitution: 

A A = 0, (64a) 
V eB [g] [AB] = ^ {AB) W B} (64b) 
TZ AB := Rf B + V e[A [W] B] - \aW A W B - \m 2 M AB = 0, (64c) 



where the density is y/^g = y —det(g AB ), the mass tensor: 

M AB = g [AB] - g C Ag B Dg [CD] + \9ba9 [CD] 9[cd\ (65) 

has been defined, and the tensor appearing in the compatibility equations is: 

^cab = -la(g A[c g D ]B + \9ab9[cd])9 {DE) W e . (66) 

One must be careful to treat totally antisymmetric derivatives properly (the structure con- 
stants now come into the curl of a vector), and translate the metric density properly. 

In order to define the action, one should note that the inverse of the metric is now: 
9ab9 BC = g CB gBA = $ A , and the direct translation of the density results in: \f—g — > 
y/ —EgE 1 = Ey/^g where g = det(g AB ) and E = det(-E^). Then ([!]) is rewritten: 

S = J m d 4 x E {-g AB R7 B - g AB V e[A [W] B] + l A A A + \ag^W A W B + \m 2 ^g [AB] ) . 

(67) 

(Note that in a Lorentz basis, the inverse of the metric is not rj.) 

Deriving the equations of motion from this action should be approached with care. As 
it stands there are too many fields (the metric and the vierbeins share degrees of freedom) 
and one typically must choose either a coordinate basis (as in Section [0]), a Lorentz basis 
(so that all symmetric metric degrees of freedom are contained in the vierbeins), or a well- 



defined combination of the two. One must also realize that (67) as it stands assumes that the 
connection T is torsion-free a priori, so that when varying the vierbein, rn must be varied 
as well. As an alternative, one may impose the torsion-free condition through additional 
Lagrange multiplier terms: L^T^q in the action, varying the full connection coefficients 
and vierbeins separately. 



V. APPROXIMATION OF THE WYMAN SECTOR SOLUTION IN A VIERBEIN 

BASIS 

In general, the spherically symmetric Killing vector analysis for a (0, 2) tensor gives 
both t — r and 9 — <fi skew components. However it is possible to show from the general 
spherically symmetric field equations that it is consistent to put either (or both) of these skew 
components to zero separately, since in either case one loses the corresponding field equation, 
and the system of equations remains consistent. Whether it is physically reasonable to do 
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this or not depends on the details of the matter coupling in the theory, and how it alters 
the global behavior of the skew sector. Here will be considered the field equations for what 
will be referred to as the Wyman sector |41]] (keeping just the 9 — cf> sector), although the 
asymptotics of the t — r sector will be discussed briefly at the end of this section, where it 
will be argued that there are no static solutions with asymptotic behavior that is dominated 
by Schwarzschild (or equivalently, Newtonian) effects. This will allow an analysis of the 
perturbation equations for the spherically symmetric modes, in order to see the effects of 
the antisymmetric background. 

In a coordinate basis, the Wyman metric looks like: 

Ifl^l = diag{7(r), -a(r), -r 2 , -r 2 sin 2 (6»)}, # [23 ] = f(r) sin(0). (68) 

(An appropriate coordinate system has been chosen in order to remove the symmetric t — r 
metric component, and fix the 9 — 9 component.) Introducing the usual choice of vierbein 
(using the functions defined by: F = f(r)/r 2 ,E Q = 1/ Jj(r), E\ = l/ya(r)): 

E/ = diag{£ ° = Eo, = E u E 2 = ±, E 3 3 = — -^}, (69) 
the metric becomes: 

\9(AB)\ = Vab, 9[23]=F, (70) 



and the density y/—g = \/T+ F 2 . 

At this point one can invert the compatibility conditions and compute the field equations 
using the method of Section [TV], given in some detail in Appendix [A|. No attempt will be 
made here to solve the field equations exactly, although numerical evidence for the existence 
of an exact solution with asymptotic behavior that matches that given here has been found 



42[, ensuring that the approximations given come from a global solution. Instead, an 
approximation will be given that describes the asymptotic behavior of the exact solution. 
The idea will be to consider the skew sector as a small correction (of order some small 
dimensionless parameter k, to be explicitly defined later) to the Schwarzschild solution far 
enough away from the source. This should be reasonable since one expects from the results of 
the perturbation in Section |T1| that the skew sector will behave asymptotically as a decaying 
exponential, while the symmetric sector should behave as ~ 1/r, so that far enough away 
from the gravitational source the skew sector should be completely dominated by GR effects. 

To lowest order in k (the skew sector) the work is already done, as the field equa- 
tion for the skew function will be essentially the same as the static perturbation about a 



Schwarzschild background already considered in Section |H4 In the vierbein basis, this is 
derived as before from 7^[ 23 ] (|A8|) , and gives: 

9 A 9 A 

Ad 2 r [F] + (A' + —)d r [F] - -(A' + -)F - m 2 F = 0, (71) 

and it is trivial to see that when using: F = f/r, this reduces to the static limit of (p2|), 
giving the asymptotic form for F: 

-mr* 

F ~ F . (72) 

mr 
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One now must consider how the presence of the skew sector affects the symmetric sector, 
particularly whether it really is a higher order effect. The asymptotic form of these correc- 
tions due to F may be calculated by considering order k 2 corrections to the vierbeins (order 
k terms will not depend on F, and so will be solely 8M S corrections), calculated from the 
symmetric field equations with F from ([72] ) acting as a source. Writing the corrections to 
the vierbeins as: _E ,i — > £q,i + Eol, an d the corrections to the field equations as TZ^ 2 \ one 
calculates: 

K { S + nf} = -—d r [VAE^ + - AFF" - -A{F'f - —FF' = 0, (73) 
r J A 2 r 



which, after translating it into a differential equation in r* and keeping only the asymptoti- 
cally dominant terms, results in: 

d^v^f + ~ 4(F ) 2 ^__. (74) 



^ J 4 

this integrates to give (the constant of integration is ignored as one could eliminate it through 



an appropriate choice of gauge as in Section III) 



771(2) r -2mr* 

+ ^ ~ l<F f—. ,75, 

Considering next 7^33 and using (|73|) leads to the asymptotic equation: 

d r *[rVAE{ 2) ] ~ -f (F ) 2 e- W . (76) 

The solution of this combined with the results of (|75|) gives (once again the constant of 
integration is ignored, this time as it would be interpretable as a perturbation of the mass 
parameter and not due to the effects of the skew sector): 



where the dominant correction to the symmetric sector is E\ , and Eq is down by o(l/r*). 
It is not hard to see that these corrections are indeed of an order higher than the effects 
in the skew sector. Clearly one may define a small parameter k = _Foexp(— mrj), where Tq 
is chosen such that F(tq) 1, to define the small size of the skew sector when r* > Tq. 
The corrections to the symmetric sector are seen to be of order k 2 , and will therefore be 



neglected in the approximation of the background required in Section [VI. 

One may attempt to do the same sort of analysis keeping the g^ij component, however 
the linearized field equation implies immediately that the field must vanish (it is identical 
to (|3"0|)). Considering higher orders in the field in an attempt to generate a solution other 
than this trivial result, the third order correction gives (writing W(r) = \jx(r)): 



/v-[01] 



Ad 2 [X] + (A' + —)d r [X] + (— + — + -m 2 )X + — + 3m 2 



0. 

(78) 
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Writing X = Y/r and transforming to the r* coordinate as before gives in canonical form: 

194 34' 3 

d 2 r , [Y] + A(— + — + ^-m 2 )Y + A{12A + 3m 2 r) = 0, (79) 
r r 2 

and keeping the dominant terms: 

d 2 r , [Y] + \m 2 Y + 3mV = 0, (80) 
easily giving the asymptotic form of the solution: 

W 2 (r*) = -2 + — cos(v/fmr*) + — sin(\/fmr*), (81) 

(where (a, 6) are arbitrary constants). The dominant part of this solution implies that W is 
imaginary, and must be discarded. However one also sees that this solution is not in fact a 
small correction to the Schwarzschild metric asymptotically, and would have to be discarded 
for that reason alone. This is not surprising as one is trying to match a function that is 
small asymptotically (by hypothesis) to one that is constant, so keeping higher orders in W 
will not change this. This implies that nontrivial static solutions that include this sector fail 
to be dominantly Schwarzschild for large r. This of course does not exclude solutions with 
asymptotic behavior that is of some other form, nor can one exclude the possibility that W 
is nonvanishing only inside some finite radius. 

VI. SPHERICALLY SYMMETRIC PERTURBATION ABOUT A WYMAN 

BACKGROUND 

In an attempt to consider the perturbation equations for NGT about a general non- 
symmetric background, one finds that the compatibility conditions prevent one from for- 
mulating the inversion in a useful form. This means that a fairly straightforward covariant 
formulation (like that given in Section |J) is n °t feasible, and instead one must treat each 
situation separately, in this case a spherically symmetric perturbation about the approxi- 
mated mNGT Wyman solution given in the previous section. Here it is demonstrated that 
despite the remaining gauge freedom in the symmetric sector, both symmetric functions 
will in general pick up time dependence from the skew sector. Although this cross cou- 
pling is demonstrated explicitly in a perturbative scenario, it will certainly persist in a more 
general sense. The results here will show that the perturbations in the symmetric sector 
pick up time dependence that is algebraically determined by the skew function F, without 
themselves becoming independent degrees of freedom. The canonical analysis of the general 
spherically symmetric system will address rigorously how many degrees of freedom exist in 
each sector. If there are more in the nonperturbative theory, one can examine the dynamical 
approach to an asymptotically flat spacetime looking for possible singular behavior similar 
to that found in 0. 

The perturbation of the Wyman metric (|70|) in a coordinate basis will look identical to 
(using the gauge choice to simplify it as before). The background vierbeins will be the 
same as those in the Wyman solution (|S9|) , where now the perturbations of the vierbeins 
and skew metric functions are related to perturbations in the coordinate basis by: 
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SW= kmM , 5F= 5k ^y\ 5E = - l ~ h ^l\ 8E 1 = - 1 - h ^4. (82) 
sja(r)-f(r) r 2 7<j) 5 2 «( r ) 5 

In the vierbein basis the metric perturbation has nonvanishing components: 

V] = SW, h m = 8F. (83) 

The approximation of the background Wyman solution given in Section [V] greatly simplifies 
the algebra necessary to develop the perturbation given in Appendix |B|. Approximating the 
symmetric sector by the Schwarzschild solution and the antisymmetric sector by (|72|), first 
order in this static antisymmetric background is kept, as is the first order in the pertur- 
bations. As one shall see, this will be a reasonable approximation since it will be possible 
to keep the perturbations small compared to the background by an appropriate choice of 
integration constants (similarly to 5M S in the Schwarzschild case). 

The field equation: ^oi] = 0, yields precisely the same field equation as in the 
Schwarzschild case ([30]), allowing one to immediately set 5W = 0. The symmetric part: 

(01) 



1 7l(m\ = 0, can be written total time derivative: 



^(oi) = d t 



2 8E 1 + %F'8F + (- - ^t)F5F + Fd r [5F] 



0. (84) 



rVA 2 V 2A' 

This last field equation is then integrated, introducing an arbitrary static function SE(r): 

2 (5E 1 -5E) = --F'5F-(--^-)F5F-Fd r [5F]. (85) 



rv 7 ^ ' 2 V 2A 

Now computing (Tt[H A b] ■— 7^oo + Tin + 27^2 ): 

Tr[*R, AB ] = ^drlrVASE^ - ^d 2 [5F] + 3AFd 2 [8F] 
r z A 

+ (—^-F + 2A'F + 5AF')d r [8F] + (3AF" + ~^-F' + 2A'F' - m 2 F)8F = 0, (86) 
and inserting (|85|) gives: 
TtFKab] = ^d r [rVA8E] 



(1 2A 2 A' 2A \ 

jd 2 [5F] - Ad 2 r [8F] - (A' + ^)d r [5F] + + + m 2 )F5F J = 0. 



(87) 



Also, 



^oo + 'Tin = — d r [VA8E + <)L> 



r si A 

F 2F 2F' 

+ -d 2 [8F] + AFd 2 [8F] + A(— + 3F')d r [8F] + A( + F")8F = 0, (88) 

A r r 

gives the equality of spatial derivatives of \ZA5E and 8Ei/ \J~A up to order F. This will be 
useful when considering: 
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[23] 



F F' 

r 2 



5 Ex 



- VA6E 



OA O A 

AF" + (A + —)F' --(A'+ -)F 

ry ry ry 



SE 1 



A 



2 A 



,2 A' 2 A 



+- -d 2 t [5F] - Ad 2 [5F] - (A' + — )d r [5F] + ( + — + m 2 )6F 



2 \A 
F F' 
r 2 



^ - VASEo 
A 



— m b - 



2A 



A 



2A' 2A 



+- -d [5F] - Ad [5F] - (A' + —)d r [5F] + (— + — + m >)5F = 0, (89) 
2 \A r r r z J 

where use has been made of fl7jj) . One derives a simple field equation by inserting (p9|) in 
(|87|) and dropping the resulting terms that are of second order in the background skew field 
F: 



d r [ry/~A5E] = 0^5E 



rVA' 



(90) 



where the constant of integration has been identified with the GR-like perturbation of the 
Schwarzschild mass parameter. 

Now fl8"8| ) can be used to replace 5E with 8E\ at this order, and Q55D to replace SEi 
with SE to find: 



[23] 



1 9 A 

-rd 2 [5F] - Ad 2 [5F] - (A + — )d r [5F] + 

-A(F )d r [-= Q) 

r J A VA 



.2 A' 2 A 



+ — + m 2 )5F 



(91) 



and using fl90"D in this yields the wave equation for SF: 

1 94 94 9 A' 

-d 2 [5F] - Ad 2 .[5F] - (A' + —)d r [8F] + (m 2 + + — )5F 

ry ry£ ry 

= — ^m 2 F + 4 - -- )(- + A . 

r 2 r 



(92) 



Note that this is a static source and so will not in itself induce any wave solutions, but as 
before the effects of a matter source will show up asymptotically The static part of the 
solution may be derived using the methods in Section 



5F 



-2F 



<5Af, 



-2mr 



*ln( 



2M, 



(93) 



and is consistent with the static solution ( [72]) derived about a Schwarzschild background 
with mass parameter M S + 5M S . Time dependent solutions are identical to those found from 
([TT|), and will induce time dependence in the symmetric sector through ([85]). Since 5E\ is 
related to 5F locally, it is not an independent degree of freedom, and since the skew field is 
short-ranged, it will not radiate energy at infinity. 

Using (|9~2"D and (|7I|) , one reduces ( |38|) to an algebraic relation for SE : 
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9 4 A F 3 A' 9 A A 

—d r [VA5E + ^=}+ AFdl[5F\ + + — F + ^F')9 r [5F] 

+K- 4 ' + |)f - + ^ + ^ + j™ 2 )^ = o- («) 

The solution for 5i?Q can be written as: 

VA5E = VA5E -^L + B{t) (95) 

where B(t) is an arbitrary function of time (removable in the usual way using the remaining 
gauge transformation), the second term corresponds to the static Schwarzschild perturbation 
from Section |HI| , and 5E solves the remainder of (|94|). Note that although not independent 
degrees of freedom, neither 5Eq nor 5E\ is static. This is in fact what one would expect when 
considering the effect of a spherically symmetric matter field to which Birkhoff's theorem 
does not apply, on the GR background. The presence of the non-static field will induce 
time dependence in the gravitational fields, without exciting any independent modes. This 
is expected to continue to be the case in NGT: the general spherically symmetric system 
should only have degrees of freedom in the skew sector. 



VII. CONCLUSIONS 

The asymptotic behavior of the antisymmetric sector for the case of a static Wyman- 
type metric has been determined, and the corrections to the symmetric sector shown to be 
negligible provided one considers regions of spacetime far enough away from the gravitational 
source. It has also been determined that if one keeps the antisymmetric t — r component, 
then one cannot have asymptotic behavior that is dominated by the Schwarzschild metric, 
and so it must be discarded. This analysis was facilitated by the introduction of a vierbein 
basis, although the formalism has been given for a general basis for completeness. 

By considering a spherically symmetric perturbation of the Schwarzschild metric, it has 
been shown that NGT does not have a rigorous Birkhoff theorem as the antisymmetric sector 
will not remain static in general. (This has also been noted previously in a Unified Field 
Theory based on Lyra geometries P3|| .) Perturbing an approximate Wyman background in a 
vierbein basis has shown that the symmetric sector is also not static in general, although no 
additional modes become excited. This is important phenomenologically since one cannot 
consider the static solutions (Schwarzschild and Wyman) as the only spherically symmetric 
exterior solutions to the field equations, and one must therefore match an interior solution 
to a non-static exterior in general. 

Perturbations of GR backgrounds have been shown to have good asymptotic behavior in 
general, since the ghost modes do not become excited and the remaining degrees of freedom 
are short ranged by construction. However this is not good enough since one expects that 
the physically interesting solutions to mNGT will not be the purely GR solutions, and one 
would therefore like to examine the behavior of perturbations on generic, asymptotically- 
flat, mNGT backgrounds. A covariant perturbative scheme, although possible in principle, 
would seem to be too complex to be of any practical value. Instead one may treat each 
case separately and consider the behavior of (perhaps several) modes about a particular 
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background, as was done here for the spherically symmetric perturbation about a Wyman 
background. 

However this also may not be adequate to fully understand the dynamics of the skew 
sector in mNGT. The lack of additional gauge invariance in the skew sector may mean that 
there are more modes in the rigorous theory that will be seen in any sort of weak field, 
perturbative analysis. To determine whether or not this is the case will require a canonical 
analysis of the full theory. Partial information may be obtained by considering the full 
set of fields in a spherically symmetric system, and looking for global information about the 
behavior of the skew modes given a general coupling to external sources. This is not likely to 
be a tractable problem in a coordinate basis, and even in a Lorentz frame the field equations 
are not expected to be particularly enlightening, due to their complication alone. However 
the canonical analysis of this system will show which fields propagate in the general case, 
and allow one to get at the dynamics of the approach to an asymptotically well behaved 
spacetime. 
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APPENDIX A: WYMAN SECTOR FIELD EQUATIONS 

Here the details of the calculation of the Wyman field equations are given, following 
the steps outlined in Section [TV]. Beginning with A A B = r\ Q\cb\ from Q53|), one finds the 
remaining components: 

A\ = -A\ = F (Al) 

The inverse of the symmetric part of the metric is the Minkowski metric: rj, whereas the full 
inverse metric is: 
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1 


1+F 2 


l+F 2 J 



(A2) 



and is necessary in order to compute the mass tensor (^). The structure constants, and 
hence the skew components of the connection T, are then easily found to be: 

T\ w] = l -C l( ? = l -E 1 d r [HE,)} 1 
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2r U 
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(A3) 



Then A from (|55| ) may be determined: 

— Ai[23] = A 3 [i2] = A 2 [3i] = Eid r [F], 

and Q from 

fil[23]=^l(^-^ r [F]), 
^3[12] = ^2[31] = |^l9 r [F]. 



(A4) 



(A5) 



Now 



and: 



is solved for A and fl53|) to calculate the connection components (T): 
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= -cot(0). 
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The field equations that remain are: 



n 00 = E 1 d r [Tl } + Tl (T° 01 + 2Tl 1 ) + 



n 22 

7£[23] 



4 1 + F 2 
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TZn — —ExdrlV^) - 2E 1 d r [Tl 1 ] v . Uiy ~ y ^- zu -v—i^ A 1 + F 2 

1 rr? 2 F 2 
TZ 33 = E 1 d r [T\ 2 } + - + I* (I* + 4r 2 ) - 2A 3 2 A 23 + — , ^ = 0, 



4 1 + F 2 



771 r 

E x d r [kU + (r^ - 4rf 12] )A 23 + 2r 22 A 3 2 - —(f + = o 



(A8) 



It can be shown that in the absence of any skew sector altogether, the Schwarzschild field 
equations are obtained, and in the absence of the mass term one has the Wyman Field 
equations O . 
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APPENDIX B: PERTURBATION EQUATIONS 



The background field will be that given in Appendix approximated by considering 
only first order contributions from the skew sector. The spherically symmetric perturbations 
about this background are given here in detail, keeping only first order in the background 
F, and setting a = 3/4 throughout. One may begin by calculating first order corrections to 
metric quantities, first the density: 



and the inverse of the full metric metric is: 



-gF5F, 
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-SW 
2F5F -5F 
5F 2F5F 



(B2) 



The tensor 5 A has remaining components: 

5A\ = SA\ = SW, 5A% = -5A\ = -5F, (B3) 
and the perturbation of the antisymmetric connection coefficients, derived from (50b): 
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(B4) 



The last equation of ( |64| ) is now solved to determine SW in terms of metric functions to 



find: 



4 25W 
SW = ~E 1 (d r [SW} + ), 

O T 

SWi = 
5W 2 



jE d t [5W], 



SW 3 = 0. 

It is then fairly straightforward to calculate the remaining SA°: 



(B5) 



^0[01] 

<5A?[oi] 



-\SW 1 , 
-\5W Q , 



-jSW a , 
\FSWa, 



(B6) 



where a e {0, 1} here and in the following. The remaining As are: 
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5A 2{a3) = -5A 3(a2) = -\F5W a , 
5A om = 2E d t [5W)-±8W 1 , 
5A 1[10] = 2E 1 d r [5W] - \8W Q , 
SA 2[2 a] = 5A 3[3a] = -\5W a , 
5A 0[23] = -E d t [SF], 
6A m = -E 1 a r [6F]-6E 1 d r [F), 
5A 3[0 2] = SA 2 [03] = E Q d t [5F], 

6A 3[12] = -5A 2[13] = E\d r [8F] + 8E x d r [F\. (B7) 



Then Sfl can be calculated: 



5Q 0[W] = E d t [5W] - \5W U 
5n m = E 1 d r [5W] - \sw , 

diZ 2 [02] = <J"3[03] = ^OWq , 

5tt 2[12] = 50 3 [i3] = \SWi, 
5O 0[23 ] =-±E d t [6F], 

1 pi rp rp x rp i 

SQ m = -~E x d r [8F] + + - -8E\d r [F], 

2 r r 2 

8tt 3 [02] = -(M 2[03] = \E Q d t [5F}, 

SQ 3[12] = -<K2 2 [i3] = \Exdr\5F] + i^^fF]. (B8) 

One is now in a position to invert the compatibility conditions ( p6| ) to solve for the 
perturbations to the connection coefficients. First the corrections to A: 

8A° 01 = -lE d t [8W], 

2 8W 
SAl, = d r [8W\), 

o T 

8Al = 8Al, = - l -E 1 { 5 ^-d r [8W]), 

8A\ X = 8A\ X = \E Q d t [8W], 
5A° 23 = -lE d t [5F], 

_ Al Ei8F 8E X F I8E1F' lEidJSF] 

SAL = — — i 1 — — 

r r 2 r 2 r 

SA 2 03 = SAl = ±E d t [SF], 

SA 2 13 = SA 3 21 = \Exd r [8F] + \F'SE X , (B9) 



and then to V: 



8T° 01 = ST], = -^d r [SE ] - SE 1 d r [\n(E )] + §-5E d r [ln(E )]., 

5r° n = sr{ = -^[sej, 

5T 2 Q3 = 5T 2 = -STL = -ST" = \E X F'8W, 
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STL = STl 
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\E 1 Fd r [5F} + \E l F'5F, 
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(BIO) 



The remaining field equations will be: ^ocb ^(oi), ^[oi]? ^-n, ^22 



% 33 and ^[23]. The 



relevant combinations will be quoted in Section [VI 
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